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ABSTRACT 

In order to improve the fit to the optical counts, the multi- wavelength source count 
model of Rowan-Robinson (2001) is improved in two ways. First, density evolution is 
incorporated to correspond to the mergers in hierarchical models of galaxy formation. 
Secondly, the evolution of dust opacity is taken into account. The luminosity and 
density evolution parameters and the characteristic opacity of the local universe are 
allowed to vary freely and the best overall fit to the B band is determined. This 
parameter set is then translated to other wavebands, showing significant improvement 
of goodness of fit in the optical, UV and near-IR bands, whilst retaining the goodness of 
fit in the far-IR and submillimetre bands of the original model. A significant fraction 
of the improvement of fit is shown to be due to the introduction of optical depth 
history. The mean value of Ay in the local universe for the best-fitting model is 0.4. 
The goodness of fit to available star formation history data is not significantly altered 
and the form of the evolution of the luminosity function shown to be consistent with 
a semi- analytic, forward evolution model. 
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1 INTRODUCTION 

The cosmic history of star formation has been extensively 
investigated in recent times. There is a general consensus 
that the star formation rate (SFR) increases with look-back 
time to peak somewhere at 1 < z < 2, then levels off, with 
a possible decline at higher redshifts. However, the issue is 
far from closed. In particular, the location of the peak and 
the behaviour of the SFR at high redshifts are still a matter 
of extensive debate. 

The nature of the star formation history (SFH) can be 
determined using a number of different observational tech- 
niques. Perhaps the most direct method is to make use of 
observations in the rest- frame near-ultraviolet (UV), as this 
part of the spectra of galaxies with no active galactic nuclei 
(AGN) are dominated by emission from short-lived, massive 
stars (e.g. Lilly et al. 1996; Madau et al. 1996; Connolly et 
al. 1997; Madau, Pozzetti & Dickinson 1998; Treyer et al. 
1998; Cowie, Songaila & Barger 1999). It is also possible to 
observe the luminosities of nebular emission lines (especially 
Hq), which are due to reprocessed ionizing radiation from 
the young, massive stars (e.g. Hogg et al. 1998; Tresse & 
Maddox 1998; Glazebrook et al. 1999; Hopkins, Connolly & 
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Szalay 2000; Teplitz et al. 2000; Pettini et al. 2001). Neither 
of these methods, however, account for the effects of dust 
extinction. 

The importance of extinction has been shown by numer- 
ous studies that show the disagreement of different SFH di- 
agnostics, which are affected by extinction in different ways 
(e.g. Mushotzky & Loewenstein 1997; Rowan- Robinson et 
al. 1997; Cram et al. 1998; Flores et al. 1999; Adelberger & 
Steidel 2000; Sullivan et al. 2001) and from the large frac- 
tion of energy in the cosmic far-infrared and submillimetre 
backgrounds (Puget et al. 1996) . Estimates of the extinction 
correction factor are around 2-10 (e.g. Meurer et al. 1997; 
Dickinson 1998; Pettini et al. 1998; Glazebrook et al. 1999; 
Steidel et al. 1999; Hopkins et al. 2000; Teplitz et al. 2000). 
It is, however, possible that some star formation could be 
totally hidden in heavily obscured objects, which would not 
be accounted for at all by a correction to the UV/optical es- 
timates (Franceschini et al. 1997; Guiderdoni et al. 1998; El- 
baz et al. 1999; but see Adelberger & Steidel 2000). Further- 
more, it has recently been shown that the far-infrared (FIR) 
to UV flux ratio breaks down for some star-forming galax- 
ies (Bell & Kennicutt 2001), ultraluminous infrared galaxies 
(ULIRGs; Meurer et al. 2000) and Lyman Break Galaxies 
(e.g. Baker et al. 2001), which makes it difficult to determine 
the SFR by simply applying an extinction correction. 
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To avoid these problems, it is necessary to look at the 
radiation re-emitted from dust heated by the young stellar 
population, which will also reveal any obscured star-forming 
population. Since thermal dust emission peaks in the rest- 
frame fax-infrared (IR), this is done with submillimetre and 
far-IR surveys (e.g. Rowan-Robinson et al. 1997; Hughes et 
al. 1998; Blain et al. 1999b). The main drawback with this 
method is the uncertainty of the contributions of cumulative 
heating of the dust by the old stellar population and AGNs 
(Blain et al. 1999c; Brusa, Comastri & Vignah 2001). 

A very different way of approaching the issue of how to 
determine the SFH is to use source count models. Broadly 
speaking, these fall into one of two categories: the so-called 
'forwards' and 'backwards' models (Lonsdale 1996). The for- 
wards method statistically follows the growth of dark matter 
haloes forwards in time by accretion and mergers (e.g. Kauff- 
mann & White 1993; Cole et al. 1994; Blain et al. 1999c; 
Granato et al. 2000; Somerville, Primack & Faber 2001) then 
assumes a semi-analytic galaxy formation paradigm and evo- 
lutionary framework. In the backwards method, by contrast, 
a local luminosity function is evolved backwards in time, as- 
suming some form of evolution, such as a change in SFR 
(e.g. Pearson & Rowan-Robinson 1996; Bales & Edmunds 
1997; Xu et al. 1998; Takeuchi et al. 1999; Pearson 2001; 
Rowan-Robinson 2001; Xu et al. 2001). 

In Rowan- Robinson (2001; hereafter RROl), a source 
count model using a parametrized approach to the SFH 
was presented. This method models the SFH as a com- 
bined power-law and exponential, the relative strengths of 
which are allowed to vary freely. This is capable of reproduc- 
ing most physically realistic forms of the SFH and allows a 
large set of possible histories to be investigated. Source count 
predictions are then compared to observations, allowing the 
range of SFH parameters to be limited. This method is sim- 
ilar in some ways to those of Guiderdoni et al. (1998) and 
Blain et al. (1999b). 

The parametrized approach of RROl provided good fits 
to data in the submillimetre and far- and mid-IR regions 
of the spectrum (10 [im < X < 1250 vim). However, the 
best-fitting model shows a deficit of faint {B > 24, K > 20) 
optical and near-IR sources compared to observations. Also, 
the number of high redshift sources are over-predicted. RROl 
proposes that these problems could be caused by density 
evolution, or by a steepening of the faint-end luminosity 
function at z > 1. 

In this work, an attempt to solve the faint optical count 
issue by allowing a variation of optical depth with redshift 
using the Calzetti & Heckman (1999) evolution model and 
by including density evolution is presented. The objective 
is to produce an effective, multi-wavelength source count 
model, based on that detailed in RROl. 

The layout of this paper is as follows. In ^ a brief 
explanation of the RROl parametrized source count model 
is presented. Then, in the optical depth history correction 
is detailed. Section ^ discusses the method of fitting to the 
optical counts and presents the best-fitting model. Section 
^ shows the predicted counts in other wavebands and the 
background spectrum, with comparisons to available data. 
In «, the nature of the fits to source count data, the nature 
of the optical depth history, the implications for the SFH 
and galaxy evolution and comparisons to a semi-analytic 



model are discussed. Finally, the conclusions are presented 



2 PARAMETRIZED SOURCE COUNT MODEL 

A backwards source count model can be constructed using 
four basic ingredients; a cosmology, a luminosity function, a 
form of evolution (luminosity and/or density) and a spectral 
energy distribution (SED). In this work, the model of RROl 
is used as a base: a more detailed description of the model 
can be found in that paper. For computational speed, an 
Einstein-de Sitter {Q,o = 1, A = 0) Universe with Ho — 100 
kms~^Mpc~^ is examined first. The results of this are used 
to limit the free parameter space of the Lambda cosmology 
(!^o = 0.3, A = 0.7). 

2.1 Luminosity Function 

The luminosity function used is determined from the IRAS 
PSCz sample (Saunders et al. 2000). This is then fitted with 
the empirical form of Saunders et al. (1990) for a single, 
composite population: 

where C*,ci and a are constants. This behaves as a power 
law for L <^ L, and as a Gaussian in log L for L ^ L^, 
L, being the 'knee' of the curve. For a and C, , fixed values 
of a = 1.09 and C. = 0.027(ffo/100)^ have been assumed. 
The other parameters take into account the change of the 
60- nm luminosity function as the rate of evolution is varied 
(see §0). 

2.2 Parametrized form of Evolution 

Changes in the SFR within a galaxy will affect its luminos- 
ity. RROl assumes that the changes of SFR are the dominant 
factor in luminosity evolution, and parametrizes this evolu- 
tion to allow a wide range of model star formation histories 
to be explored. The form of the star formation rate, cf>t{z) (in 
units of Mq yr~^ Mpc~^), adopted to allow most physically 
realistic, single population SFH scenarios to be reproduced 
is: 

ili^=fexp0(l-l)l (1)", 

where P and Q are positive variable parameters (see fig. 1 
of RROl for examples of this form). In an Einstein-de Sitter 
(E-dS) Universe, this becomes: 

= {exp Q [1 - (1 + .)-«/^] } (1 + zy-'^/' . (1) 

The physical meanings of the SFH parameters are as 
follows. P measures how steeply the process of buildup of 
rate of star formation due to mergers of fragments, which are 
incorporated into galaxies, occurs with time. Q describes an 
exponential decay in the star formation rate, on a timescale 
to/Q, which is a product of the process of exhaustion of 
available gas for star formation - this being the result of 
a competition between the processes of gas being formed 
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into stars and it being returned to the ISM by winds and 
supernovae. 

It is assumed, for computational purposes, that </>*(2) = 
for 2 > 10. Although this introduces an unphysical discon- 
tinuity in eqn . (2.2), RRO 1 finds that the effect on the results 
is negligible. 

Now that the luminosity evolution has been selected, 
the remaining parameters in the luminosity function can be 
approximated in the E-dS case by least-squares fitting the 
Saunders et al. (1990) form to the Saunders et al. (2000) 
data set as follows: 

a = 0.744 - OmSW, 

IgL, = 8.50 + 0.05W- 2 Ig (Ho/100), 

where 

W = Q - 1.125P. 

In a hierarchical galaxy-formation scenario, it is ex- 
pected that a lower total density of galaxies will be seen 
at low redshifts as a result of merger events. In addition to 
the luminosity evolution, this density evolution also needs 
to be considered. This can be parametrized simply as: 

p(2)=p(0)(l + 2)". (2) 

In this case, the evolution parameter, W, becomes 
W = Q - 1.125(P-0.57i). 

Many source count models consider only pure luminos- 
ity evolution or pure density evolution, but here both forms 
will be considered simultaneously. 

In the Lambda cosmology with density evolution, these 
parameters become: 

a = 0.783 - O.OSOVy, 

Ig L* = 8.39 + 0.05 W - 2 Ig {Ho /lOO) , 

where 

W = Q - 1.067(P--0.5n). 

2.3 Spectral Energy Distributions 

To transform the 60- (xm luminosity function to other wave- 
lengths, some model of the SED must be constructed. RROl 
explored a variety of assumptions about the SEDs. The 
method that was most consistent with available data was 
to derive the counts separately for each of four components 
- starburst (M82-like), 'cirrus' (normal, non-active galaxy), 
ULIRG (Arp 220-like) and AGN - then sum them. 

For the starburst and cirrus components, predictions for 
infrared SEDs from Efstathiou, Rowan-Robinson & Sieben- 
morgen (2000) are used, and near-IR/optical/UV SEDs are 
added, corresponding to an Sab galaxy (Yoshii & Takahara 
1988) for the cirrus component and an H ll galaxy, Tol 1924 - 
416 (Calzetti & Kinney 1992), for the starburst component. 
These are selected as 'typical' galaxies, rather than as aver- 
ages. For the cirrus component, the optical SED has been 
divided into a contribution of young, high-mass stars (dom- 
inating at A < 0.3 (xm) and of old, low-mass stars (dominat- 
ing at A > 0.3 [im). This is obviously a quite approximate 
treatment, but produces reasonable counts in the K- and 
B-bands. For the AGN component, the dust torus model of 



Rowan- Robinson (1995) is used (see also Rowan- Robinson & 
Efstathiou 1993 and Efstathiou & Rowan- Robinson 1995). 
The complete SEDs are shown in fig. 3 of RROl. 

The proportions of the four components at 60- vim as 
a function of luminosity were chosen to give correct rela- 
tions in colour-colour diagrams, as detailed in RROl. The 
relative proportion of the 60- (xm emission due to AGNs is 
derived from the 12-vi-m luminosity function of Rush, Malkan 
& Spinoglio (1993). 



3 OPTICAL DEPTH HISTORY 

The models of RROl fit the available source count data well 
from mid-infrared to submillimetre wavelengths. However, 
in the if-band, the number of faint {K > 20) sources is 
under-predicted (see fig. 18 of RROl; also plotted on Fig. ^ 
of this work) by the preferred P = 1.2, Q — 5.4 model. This 
problem is even more pronounced in the B-band, where the 
faint {B > 24) counts are significantly under- predicted (see 
fig. 19 of RROl; also plotted on Fig. hi of this work) by the 
same model. 

This faint blue excess could be caused by a number of 
factors, such as uncertainties in the faint end of the local 
luminosity function, some form of evolution or a missing 
population of objects (see review by Ellis 1997). One factor 
that is known to occur, and could cause a faint blue excess, is 
the variation of optical depth with redshift (Calzetti 1999; 
Pei, Fall & Hauser 1999). It is probable that the opacity 
changes difi'erently for the various types of sources, but we 
have used a 'characteristic' opacity variation, in the manner 
of Calzetti 1999 and Pei et al. 1999, as this reduces the 
number of free parameters. 

Star formation produces stars with a range of masses, 
some of which will be sufficiently massive to explode as su- 
pernova very quickly (~ lO" years; Leitherer & Heckman 
1995), injecting heavy metals into the interstellar medium 
and producing large quantities of dust on short time scales 
(~ 100 - 200 Myr, e.g. Jones et al. 1994; Dwek 1998). The 
global star formation rate varies with redshift, so the co- 
moving density of metals, flm, and hence the dust density, 
will also vary. This variation takes the form of an increase 
with redshift, peaking somewhere between z ~ 1 and z ~ 3, 
then declining (see fig. 8, bottom panel of Pei et al. 1999). 
This form is the result of a competition between the release 
of metals into the interstellar medium (ISM) and the metals 
being taken-up from the ISM as the gas is locked into stars. 

As the opacity of the ISM is determined by the com- 
bination of the dust column density (dependent on Qm and 
the gas column density), the extinction curve and the dust 
distribution, it is possible to model the opacity change with 
redshift. This new factor can be included in the source count 
model of RROl by altering the luminosity distance accord- 
ingly. 

Initially, the obscuration model of Pei & Fall (1995) 
was investigated, using the dust model of Rowan- Robinson 
(1992) and the colour excess model of Pei et al. (1999), 
assuming the local mean fraction of starlight absorbed by 
dust, Avq ~ 0.4. Here, the obscuration of the flux, F, is 
parametrized as: 



Fobs{\) = Fo{\)W 



(-0.4x0.4(Es(B-V))k{\)/3.1) 
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where Es{B ~ V) is the colour excess and fc(A) is the dust 
obscuration factor. 

However, this model does not solve the faint excess 
problem, as the opacity drops at 2; > 1, and we are trying 
to reduce the number of faint (hence mostly high redshift) 
galaxies. A model is required which predicts that the optical 
depth does not drop until a higher redshift. In order to meet 
this criterion, the Calzetti (1999) obscuration model was in- 
vestigated, using model 3 from Calzetti & Ifeckman (1999) 
for E{B — V), as this does not begin to drop until z > 3 (see 
fig. 9, top panel of Calzetti & Heckman 1999). The Calzetti 
(1999) model analytically parametrizes the 'net' obscuration 



i^.fc.(A) = Fo(A)10-«-^^-(^-^"=(^\ 

based on observations of high redshift galaxies, where: 

k{X) = 2.656 ( -2.310 + j + 4.88 

for 0.63 [im ^ A 1.60 (xm, and 



(3) 



fc(A) = 2.656 -2.156 



1.509 
A 



0.198 0.011 \ 



for 0.12 ^.m s; A < 0.63 ^.m. 

Outside this range of A, the Rowan-Robinson (1992) 
dust form is used. This is based on a multigrain model, fit- 
ted to the interstellar extinction curve and the far-IR spec- 
trum of dust in our galaxy. A look-up table of values of 
A\/E[B — V) for different wavelengths is used, interpolat- 
ing where necessary. 

This correction for the optical depth history requires a 
new parameter, Avq, to be introduced into the source count 
model. This represents the mean V band extinction in the 
local {z = 0) universe and can be varied whilst retaining the 
dependence on redshift of the Calzetti & Heckman (1999) 
model 3 extinction, g{z). It is introduced into the net ob- 
scuration (eqn. ^ as follows: 

i^o6.(A) =Fo(A)10-"-"^^^(^', 
where 

Ax{z) ^Avo X (^^) xg{z). 

The effect of the optical depth variation introduced 
here is wavelength-dependent. The magnitude of the effect 
is smaller at higher wavelengths. It was found that at wave- 
lengths of approximately 5-|i.m and higher the effect is neg- 
ligible. For computational speed, therefore, the correction 
was applied only at A < 10 [im. 



4 FITTING OPTICAL COUNTS 

The source count model, with the optical depth history cor- 
rection, has four free parameters: the characteristic visual 
extinction in the local universe, Avg, the density evolution 
parameter, n, and the two star formation parameters, P 
and Q. To determine if the opacity correction could solve 
the fitting problem in the B band, these parameters were 
varied and the resulting model predictions compared to the 
available observations. 

For a given waveband, there are two types of data avail- 
able against which the models can be compared: source 



count data and redshift distribution data. Generally, it is 
possible to obtain large statistics for source counts much 
more easily than for redshift distributions (few spectro- 
scopic redshifts are available and photometric redshifts can 
be unreliable). As the source counts are much more robust, 
statistically speaking, the procedure adopted was to first 
use the available redshift data - a catalogue of photomet- 
ric and spectroscopic redshifts in the Hubble Deep Field 
(Fernandez-Soto, Lanzetta & Yahil 1999) derived from data 
at 0.3, 0.45, 0.6 and 0.8 [im (filters F300W, F450W, F606W 
and F814W, respectively) and in the J, H and TiT bands (1.2, 
1.6 and 2.2 |i.m, respectively) - to constrain the parameter 
space to those sets that were as good as, or provided an 
improvement over, the P = 1.2, Q = 5.4 E-dS model of 
RROl. This was done by generating a large number of mod- 
els (~ 2500) with a wide range of different parameters, then 
using a reduced method to remove those that have a 
worse fit than the P = 1.2, Q = 5.4 model. 

As the number of sources at high redshift for the avail- 
able data is small, the bins z — 4.0 - 4.5, 4.5 - 5.0 and 
z > 5.0 are summed up to give one bin, z > 4.0. This gives 
the number of model points, rim = 12 and the number of 
free parameters, ric — i {P,Q,n and Avq) and, hence, the 
number of degrees of freedom, u — 8. Around 200 models 
had better fits than the P = 1.2, Q = 5.4 model. 

After a reasonable fit to the redshift distributions had 
been obtained, the source count data was investigated to 
find the best-fitting model. The data used are compiled from 
Anglo- Australian Telescope (AAT) observations to a depth 
of i? ~ 23.5 (Jones et al. 1991), Hubble Space Telescope 
observations of the Hubble Deep Fields (North and South) 
and William Herschel telescope (WHT) observations to a 
depth of B ~ 28.5 (Metcalfe et al. 2001), and observations 
made during the commissioning phase of the Sloan Digital 
Sky Survey to a depth of B ~ 20 (Yasuda et al. 2001). 

For the models that showed an improved fit to the red- 
shift distribution, the Xv method was again used to find 
those that had a better fit to the source count data than 
the P — 1.2, Q = 5.4 model. In this case, 4 model parame- 
ters are compared to the 28 data points, giving v — 24. All 
but six of the parameter sets examined showed a better fit 
than the P = 1.2, Q = 5.4 model. Of these, the best-fitting 
model was P = 1.5, Q = 5.4, n = 1.3, ylvo = 0.4. As can be 
seen from Figs. |l| and ^ this model gives a reasonable red- 
shift distribution whilst significantly improving the fit to the 
source count data. Since the luminosity and density evolu- 
tion parameters have been altered, these figures also show a 
model with the same evolutionary parameters, but no opti- 
cal depth correction (P = 1.5, Q = 5.4, n = 1.3, Av^ = 0.0). 
This shows that the improvement of goodness of fit is not 
due solely to the change in evolutionary parameters, but 
relies on the new optical depth history. 

Once this had been done for the E-dS cosmology, the 
Lambda model was investigated. The best-fitting Lambda 
model of RROl had SFH parameters P = 3.0, Q = 9.0. As 
the Lambda model with density evolution is very computa- 
tionally intensive, a 'first guess' was made, based on the re- 
sults of fitting the E-dS model. In the E-dS case, the SFH pa- 
rameter Q did not change and P was slightly increased when 
density evolution was introduced, so Q was fixed and P al- 
lowed to vary. As the best-fitting value of Avq was consistent 
with local observations, this was also fixed, so allowing only 
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Figure 1. Differential source counts at 0.44-|am. Data from Jones 
et al. (1991) [squares], Metcalfe et al. (2001) [triangles] and Ya- 
suda et al. (2001) [diamonds]. Models are the best-fitting E-dS 
from RROl, P = 1.2, Q = 5.4 [short-dashed curve], the best- 
fitting E-dS from this work, P = 1.5, Q = 5.4, n = 1.3, ^vb = 0-4 
[solid curve] and the comparative model P = 1.5, Q = 5.4, n = 
1.3, = 0.0 [dotted-dashed curve]. Lower curves are predic- 
tions for quasars in the P = 1.2, Q = 5.4 [long-dashed curve] 
and P = 1.5, Q = 5.4, n = 1.3, Ay^ = 0.0 [dotted-dashed curve]. 
Quasar data from Boyle, Shanks & Peterson (1988) [inverted tri- 
angles]. The goodness of fit is improved in the new model, due to 
both the changes in evolutionary parameters, and the introduc- 
tion of the optical depth correction (see text for discussion). 

the density evolution parameter, n, and P to vary freely. The 
best-fitting model was P = 3.4, Q = 9.0, n = 1.1, Avq = 0.4. 
As can be seen from Fig. ^ this model is a very good fit 
to the B band source counts, so no further variations of the 
parameters were investigated. 



5 PREDICTED COUNTS 

Once the best-fitting sets of parameters in the B band had 
been determined (P = 1.5, Q = 5.4, n — 1.3, Avq ~ 0.4 for 
E-dS, P = 3.4:, Q = 9.0, n = 1.1, Avo = 0.4 for A), these 
sets were used to generate models in other wavebands to 
determine if the fit in the UV/optical/near-IR wavebands 
had been improved. 

The best-fitting models of RROl under-predict the faint 
counts in the UV/optical/near-IR regions. If the new opti- 
cal depth history correction works, it should increase the 
goodness of fit in these other bands without tuning the free 
parameters. 

The fits to the U band is shown in Fig. ^. The overall 
goodness of fit has been improved significantly across the 
whole fiux range for both cosmologies, although the faint 
{U > 24) counts are still under-predicted in the E-dS case. 
Figs. ^, y and ^ show the fits to the R, H and K bands, 
respectively. In each case the goodness of fit has been sub- 
stantially improved for both cosmologies, though the very 
faintest counts tend to be under-predicted in the E-dS case. 




z 



Figure 2. Predicted redshift distribution at 0.44- |am, flux range 
-5.2 > lgS{Jy) > -6.8 (corresponding to 22 < B < 26). Models 
are the best-fitting from RROl, P = 1.2, Q = 5.4 [dashed curve], 
the best-fitting from this work, P = 1.5, Q = 5.4, n = 1.3, Avq = 
0.4 [solid curve] and the comparative model P = 1.5, Q = 5.4, n = 
1.3, Avq = 0.0 [dash-triple dot curve]. Data from Fernandez-Soto 
et al. (1999). Errors arc Poisson distribution values. The bin 
shown as 2 = 5.0 — 5.5 is z > 5.0. 

The Av correction is wavelength-dependent and is neg- 
ligible for A > 5 iim, so this will have no effect on source 
counts at mid-IR, far-IR and submillimetre wavelengths. 
However, the values of P and n are different to the pre- 
ferred models of RROl, which will introduce changes. Pre- 
dictions at 15, 60, 170 and 850 \im are shown in Figs. ^, ^ 
[l^ and |l^, respectively. Across all these wavebands, and for 
both cosmologies, the best-fitting models of this work pre- 
dict higher values of the faint source counts than the RROl 
models. For 60 and 170 [im, the available observations are 
not deep enough to be able to distinguish between the mod- 
els. At 15 and 850 \im there is no significant change in the 
overall goodness of fit. 



6 DISCUSSION 
6.1 Fits to Data 

The P = 1.5, Q = 5.4, n = 1.3, Avq = 0.4 model was se- 
lected to provide the best possible fit to the B band source 
counts (Fig. |l|), whilst not worsening the fit of the redshift 
distribution (Fig. As expected, this model provides a sig- 
nificant improvement of fit in this band. It does, however, 
still slightly under-predict the faintest end {B > 26) of the 
source count slope. Models with a higher value of Avq can 
provide a better fit to the faint end, but significantly de- 
crease the goodness of fit to the redshift distribution. The 
P = 3.4, Q = 9.0, n = 1.1, Avg = 0.4 model in the A cosmol- 
ogy, however, fits the faint end of the counts well. 

Having introduced the optical depth variation correc- 
tion to make the model more physical, and selected the 
parameters to fit the B band, the model is then 'trans- 
lated' to other wavebands. The fits in the optical, UV 
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Figure 3. Differential source counts at 0.44- mm in a A cosmology. 
Data from Jones et al. (1991) [squares], Metcalfe et al. (2001) 
[triangles] and Yasuda et al. (2001) [diamonds]. Models are the 
best-fitting from RROl, P = 3.0, Q = 9.0 [dashed curve] and the 
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Figure 5. Differential source counts at 0.64- mm. Data from Jones 
et al. (1991) [diamonds], Metcalfe et al. (1991, 2001) [squares, 
triangles]. Models are the best-fitting E-dS from RROl, P = 
1.2, Q = 5.4 [dashed curve] and the best- fitting E-dS from this 
work, P = 1.5, Q = 5.4, n = 1.3, Ay^ = 0.4 [solid curve]. 
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Figure 4. Differential source counts at 0.36- ^m. Data from Koo 
(1986) [squares], Jones et al. (1991) [diamonds] and Metcalfe et 
al. (2001) [triangles]. Models are the best-fitting E-dS from RROl, 
P = 1.2, Q = 5.4 [dashed curve], the best-fitting E-dS from this 
work, P = 1.5, Q = 5.4, n = 1.3, Ay^ = 0.4 [solid curve] and the 
best-fitting A from this work, P = 3.4, Q = 9.0, n = 1.1,j4vq = 
0.4 [dotted-dashed curve]. 



and near- IR bands ( U, R, H, K; Figs. |, |, |, |, respec- 
tively) are all substantially better for both cosmologies, al- 
though there is again an under-prediction in the faint ends 
{U > 24, 7? > 25, H >21,K > 21) of the slopes for the E-dS 
model. 

Quasar counts in the optical {B band) region are shown 
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Figure 6. Differential source counts at 1.6-|am. Data from Met- 
calfe et al. (in preparation) [diamonds] and Yan et al. (1998) 
[squares]. Models are the best-fitting E-dS from RROl, P = 
1.2, Q = 5.4 [dashed curve] and the best- fitting E-dS from this 
work, P = 1.5, Q = 5.4, n = 1.3, Ay^ = 0.4 [solid curve]. 



in Fig. There are, however, insufRcient faint counts to be 
able to distinguish between the models. 

As the parameter set used is derived from optical ob- 
servations, it is possible that the fit in the mid-IR to sub- 
millimetre bands might be worsened. Since the goodness of 
these fits is the main strength of the RROl model it is im- 
portant to check the effects of altering the parameter set in 
these wavebands. At 60 and 170 \Lm (Figs. |and0), insuf- 
ficient faint source counts are available to be able to distin- 
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Figure 7. Differential source counts at 2.2- (om. Data from 
Mobasher, Ellis & Sharpies (1986) [squares], Gardner et al. (1996) 
[diamonds], Huang et al. (1997) [crosses] and McCracken et al. 
(2000) [triangles]. Models are the best-fitting E-dS from RROl, 
P = 1.2, Q = 5.4 [dashed curve], the best-fitting E-dS from this 
work, P = 1.5, Q = 5.4, n = 1.3, = 0.4 [solid curve] and the 
best-fitting A from this work, P = 3.4, Q = 9.0, n = 1.1, = 
0.4 [dotted-dashed curve]. 



Figure 9. Differential source counts at GO-yim. Data from Lons- 
dale et al. (1990) [triangles]. Rowan- Robinson et al. (1990) [di- 
amonds] and Gregorich et al. (1995) [squares]. Models are the 
best-fitting E-dS from RROl, P = 1.2, Q = 5.4 [dashed curve], 
the best-fitting E-dS from this work, P = 1.5, Q = 5.4, n = 
1.3, Avq = 0.4 [solid curve] and the best-fitting A from this work, 
P = 3.4, Q = 9.0, ra = 1.1, Ay^ = 0.4 [dotted-dashed curve]. 
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Figure 8. Differential source counts at 15- [am. Data from Altieri 
et al. (1999) [diamonds], Elbaz et al. (1999) [squares] and Serjeant 
et al. (2000) [triangles]. Models are the best-fitting E-dS from 
RROl, P = 1.2, Q = 5.4 [dashed curve], the best-fitting E-dS 
from this work, P = 1.5, Q = 5.4, n = 1.3, Ay^ = 0.4 [solid curve] 
and the best-fitting A from this work, P = 3.4, Q = 9.0, n = 
1.1,j4vo = 0.4 [dotted-dashed curve]. 
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Figure 10. Integral source counts at 170- |jm. Data from Dole 
et al. (2001). Models are the best-fitting E-dS from RROl, P = 
1.2, Q = 5.4 [dashed curve], the best-fitting E-dS from this work, 
P = 1.5, Q = 5.4, n = 1.3, Ayg = 0.4 [solid curve] and the best- 
fitting A from this work, P = 3.4, Q = 9.0, n = 1.1, ^vb = 0.4 
[dotted-dashed curve]. 

guish between the models. At 15-\im (Fig. H) the goodness 
of fit is improved in the faint (S < 1 mjyj region for the 
E-dS cosmology and essentially unchanged for the A case. 
At 850- ^tm (Fig. |l^) the goodness of fit is decreased in both 
cosmologies, but the models are still within the limits of the 
observations. 
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Figure 11. Integral source counts at 850- |jm. Data from Hughes 
et al. (1998) [cross], Blain et al. (1999a) [squares], Eales et al. 
(1999) [inverted triangle]. Fox (2000) [triangle] and Scott ct al. 
(2002) [diamonds]. Models are the best-fitting E-dS from RROl, 
P = 1.2, Q = 5.4 [dashed curve], the best-fitting E-dS from this 
work, P = 1.5, Q = 5.4, n = 1.3, Ay^ = 0.4 [solid curve] and the 



best-fitting A from this work, P 
0.4 [dotted-dashed curve]. 



3.4,Q = 9.0, n = 1.1, A 



The best-fitting model was also compared to the avail- 
able data on the integrated background spectrum (Fig. 
Here, the data is compiled from surveys in the optical and 
UV (Pozzetti et al. 1998), mid-IR (Serjeant et al. 2000) 
and far-IR and submillimetre (Puget et al. 1996; Fixsen et 
al. 1998; Lagache et al. 1999). The new E-dS model [solid 
curve] shows a better fit in the optical, UV and near-IR re- 
gions, but a worse fit in the mid-IR (15- vim), although this 
point is derived by summing counts. The goodness of fit in 
the far-IR (140, 240, 350 and 500 |i-m) is approximately the 
same, but the submillimetre (750 and 850 ^.m) background 
is over-predicted. The dotted curve in this figure shows the 
comparative P — 1.5, Q = 5.4, n — 1.3, Avg = 0.0 model 
to illustrate that it is the optical depth correction, not just 
the change in P and n that is responsible for the better fit. 
The A model is little different to the best-fitting A model of 
RROl - the submillimetre and FIR are over-predicted, but 
close to the upper limits of the errors. 

The over-prediction of faint submillimetre source counts 
is dominated by the introduction of density evolution. There 
are certain limitations to the model that might, if corrected 
for, give a good fit to the optical region with a lower density 
evolution parameter. Possible candidates for improvement 



(i) Accounting for ellipticals. The count model is based on 
the 60-ii.m luminosity function, which provides a good fit in 
the infrared. However, this does miss out elliptical galaxies, 
which have little or no 60-(i.m emission. A simplistic way of 
remedying this would be to multiply the old star component 
by a factor to account for the missing ellipticals (see §|6.2|). 

(ii) Separately evolving SED components. In the current 
model, the form of evolution adopted is applied to all the 



g A (/xm) 

Figure 12. Predicted spectrum of integrated background. Data 
from Puget et al. (1996) [crosses], Fixsen et al. (1998) [diamonds], 
Pozzetti et al. (1998) [triangles], Lagache et al. (1999) [squares] 
and Serjeant et al. (2000) [inverted triangle]. Models are the 
best-fitting E-dS from RROl, P = 1.2, Q = 5.4 [dashed curve], 
the best-fitting E-dS from this work, P = 1.5, Q = 5.4, n = 
1.3, Avo = 0.4 [solid curve], P = 1.5, Q = 5.4, n = 1.3, = 0.0 
[dotted curve] and the best-fitting A from this work, P = 3.4, Q = 
9.0, n = 1.1, Ayjj = 0.4 [dotted-dashed curve]. Note that the fit in 
the optical is much improved, but the submillimetre counts are 
over-predicted (see text for discussion of ways to resolve this) . 



components equally. This simplifies the model and produces 
reasonable fits to source count data, but is not ideal. Differ- 
ent types of galaxies will evolve differently with redshift, so 
separating them out would make the model more physical. 

(iii) Redshift-dependent cirrus SED. The RROl model 
does not account for the fact that the temperature of dust, 
and hence the SED, of cirrus galaxies will evolve with red- 
shift. Including this, perhaps as a look-up table of SEDs at 
different redshifts, would again make the model more phys- 
ically realistic. 

Introducing one or more of these might allow the density 
evolution parameter, n, to be lowered, and hence the fit to 
the faint submillimetre counts improved, whilst retaining 
the goodness of fit in the optical region. The first of t hese 
factors is quite simple to correct for and is discussed in §6.2, 



6.2 Ellipticals 

Initial results from the IRAS sky survey showed that very 
few elliptical galaxies had been detected in the FIR (de Jong 
et al. 1984). Further 

Analysis of the IRAS sky survey showed that ^ 35 per 
cent of ellipticals are detected at 60- vim to a ffux densi- 
ties limit of ~ 300 mjy (Knapp et al. 1989). As the 60-ii-m 
emission is dominated by emission from cold (T ~ 20-40 
K) dust, this implies that ellipticals typically have low dust 
densities (e.g. Goudfrooij & de Jong 1995; Merluzzi 1998). 
These galaxies will not be accounted for by the 60-vim lumi- 
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Figure 13. Differential source counts at 0.44-mm. Data from 
Jones et al. (1991) [squares], Metcalfe et al. (2001) [triangles] 
and Yasuda et al. (2001) [diamonds]. Models are the best-fitting 
from this work, P = 1.5, Q = 5.4, n = 1.3, Ay;, = 0.4 before 
[solid curve] and after [dashed curve] modification to account for 
ellipticals. The effect of accounting for ellipticals is negligible. 



nosity function, which is used by the RROl model, so some 
correction is needed. 

Martel, Premadi and Matzner (1998) have used the ob- 
served galaxy morphology-density relation (Dressier 1980) 
and a Monte Carlo technique based on the distribution of 
dark matter to simulate ~ 40,000 galaxies in the local uni- 
verse. For an E-dS Universe, they find that ellipticals make 
up, on average, ~ 15 per cent of galaxies. Hence, since only 
~ 35 per cent of ellipticals are detected at 60-[im, the RROl 
model misses out ~ 10 per cent of 'cirrus' galaxies. 

A simple way to correct for this is to multiply the num- 
ber density of the low-mass cirrus component by a factor of 
1.104 at each flux bin for those models that split the cirrus 
component into low- and high- mass (A < 20ujn). The result 
of applying this correction is shown in Fig. hsl The only ef- 
fect is a very slight increase in the brightest counts, where 
the low-mass component is dominant. The overall effect of 
correcting for ellipticals is negligible. 



6.3 Optical Depth 

In order to determine how much of the improvement of fit 
is due to introducing the optical depth correction, and how 
much is due to the change of P and n, Figs. |l] and |^ show 
a model with the same parameter set as the best-fitting E- 
dS from this work {P = 1.5, Q = 5.4, n = 1-3), but with 
Avq set to zero, so that the optical depth correction is not 
applied. Fig. ^ shows that the alteration in P and n results 
in a general increase in the number of counts, especially in 
the fainter regions, whereas the addition of the change in 
optical depth is responsible for increasing the slope at the 
faint end of the curves. This shows that it is necessary to 
include the optical depth correction to achieve this goodness 
of fit. 



The overall goodness of fit to the redshift distribution 
(Fig. 1^ is little changed between the models. The change 
in shape is a result of both the optical depth correction and 
the change in evolutionary parameters. 

Observations of local spiral galaxies (e.g. Berlind et. 
al. 1997; Gonzalez et al. 1998) show the visual extinction, 
Av ~ 0.1 - 0.2 in the inter-arm regions and Av 1 in 
the arms and centres (Giovanelli et al. 1995; Moriondo, Gio- 
vanelli & Haynes 1998). The best-fitting model of this work 
has Avq ~ 0.4, which is perfectly reasonable, given known 
observational results. It is a validation of the optical depth 
history correction that the value of Avq for the best-fitting 
model is a physically realistic value purely by fitting to 
source count data, with no tuning to achieve this result. 



6.4 Star Formation History 

The main strength of the parametrized approach to the SFH 
as presented in §2.2 (eqn. 2.2), is that it is capable of repro- 
ducing most physically realistic, single-population scenarios. 
This allows a wide range of models to be explored. 

Once the best-fitting model in the B band has been 



determined, the luminosity and density evolutions (eqns. 
& H) can be combined to give a star formation history: 



2.2 



exp 



(1 



The predicted SFH for the E-dS cosmology is shown in Fig. 
compared to the P = 1.2, Q = 5.4 model of RROl. There 
is no significant change in the goodness of fit. The SFR at 
high redshifts {z > 2) is probably too high, due to the intro- 
duction of density evolution. As with the over-prediction of 
the faint submillimetre counts, this might be resolved by al- 
lowing the components to evolve separately or evolving the 
SEDs with redshift. For the A case, shown in Fig. ha, the 
goodness of fit is improved in the high-redshift region. 

It is important that the form of the star formation his- 
tory in the model produces a total mass of stars consis- 
tent with observed values. Lanzetta, Yahil & Fernandez- 
Soto (1996) used optical and UV observations of the Hub- 
ble Deep Field to determine the total mass fraction of stars, 
SI, = 0.0030±0.0009/i"^ . As this is not corrected for extinc- 
tion, this is the lower limit on the actual mass of stars. The 
upper limit is the mass fraction of baryons, which is found 
by Bania, Rood & Balser (2002) to be ^i, = 0.020t° °°l h'^. 

The theoretical value of fl* can be calculated from 

= ^ = 10-"**ft-^to<^.(to)C, 



pc 



where 



(t>*{to) 



d(t/tQ) 



The SFR in the local universe, (j>*{to), is derived from the 
60-p.m luminosity function by Rowan- Robinson (2002) using 
eqn. (7) of Rowan- Robinson et al. (1997), modified to take 
into account the latest Bruzual & Chariot galaxy evolution 
models (Madau et al. 1998), assuming that the fraction of 
optical-UV light being re-radiated in the IR, e — 2/3: 
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Figure 14. Star formation history in an E-dS cosmology. Mod- 
els are the best-fitting from RROl, P = 1.2, Q = 5.4 [dashed 
curve], and the best-fitting from this work, P = 1.5, Q = 5.4, n = 
1.3, Ayg = 0.4 [solid curve]. Data points derived from IR (Rowan- 
Robinson et al. 1997 [open triangles]; Flores et al. 1999 [open dia- 
monds]; Mann et al. 2002 [open squares]), submillimetre (Hughes 
et al. 1998 [inverted open triangle]), extinction-corrected UV 
(Gallego et al. 1995 [filled diamond]; Gronwall 1999 [filled tri- 
angle]; Steidel et al. 1999 [inverted filled triangles]; Sullivan et al. 
2001 [cross]) and re-analysed photometric redshifts of HDF galax- 
ies, corrected for extinction (Rowan-Robinson 2002, in prepara- 
tion [filled squares]). There is no significant change in the good- 
ness of fit. 



0.0 



o -0.5- 



e 



-1.0 - 



-2.0 - 







• — 











'-^ 1 


1 



^ SFH 
^Lambda) 



4 



Figure 15. Star formation history in a A cosmology. Models are 
the best-fitting from RROl, P = 3.0, Q = 9.0 [dashed curve], and 
the best-fitting from this work, P = 3.4, Q = 9.0, n = 1.1, Avq = 
0.4 [solid curve]. Data points are derived from infrared data, ul- 
traviolet data with correction for dust and photometric redshifts 
of HDF galaxies. References are the same as Fig [l4[ 



r 10 

where the 60-|i.m luminosity density, Zeo = h Jj^qs ri{L)dL. 
Hence, given to = 2/3i?o, 

f7* = 10~"'^''e^6o/i"^ (4) 



For the P = 1.2, Q = 5.4 model of RROl, ko = 
4.38 X 10^ /iLgMpc-^ and C = 5.74, giving = 0.0013 h'^. 
For the best-fitting model of this work, l^o ~ 4.58 x 
10'' /iLqMpc"^ and ^ = 12.45, giving fl, = 0.0029 h'^ . 
In this scenario, the P — 1.5, Q = 5.4, n = 1.3, Avg = 0.4 
model is most consistent with observations. 

However, the age of the Universe for this calculation has 
been set to to ~ 6.5 Gyr, which, although consistent with 
the cosmology, is not consistent with the observationally- 
determined age. Using a more realistic value of to = 13 Gyr, 
eqn. H becomes: 



f7* 



10" 



(5) 



This gives — 0.0026 h for the best-fitting model of 
RROl and n, = 0.0058 h'^ for the best-fitting model of 
this work. In this case, the RROl model is more consistent 
with the observed lower limit value of fl,. The higher fl, 
in the best-fitting model is due, as with the over-prediction 
of the submillimetre background, to the inclusion of density 
evolution. 



6.5 Comparison to semi-analytic model 

It is possible to compare the evolution from this 
parametrized backwards model against a semi-analytic for- 
wards model to determine if they are compatible. Fig. |l6| 
shows the differential luminosity function at different red- 
shifts. The semi-analytic model used is based on work by 
Kaviani, Haehnelt & Kauffmann (2002). 

The two luminosity functions are broadly similar, fit- 
ting best in the intermediate redshift range 2 ~ 1 - 2. The 
results are broadly consistent, which is interesting, given 
that they are derived using entirely different methods. A 
more extensive comparison to the semi-analytic model will 
be conducted in a future work. 



7 CONCLUSIONS 

An attempt to improve the optical count fit of the RROl 
IR-defined source count model by introducing density evo- 
lution and a variation of optical depth with redshift has been 
presented. 

(i) By introducing an optical depth correction and tun- 
ing the luminosity and density evolution parameters, the 
goodness of fit to the source counts in the B band has been 
significantly increased. This improvement is not due solely 
to the change in existing parameters, but also substantially 
affected by the introduction of the optical depth correction 
and density evolution. 

(ii) The parameter set derived in the B band can then 
be translated to other wavebands, showing significantly im- 
proved goodness of fit in the optical, UV and near-IR as well 
as reasonable fits in the far-IR and submillimetre. 

(iii) In the Einstein-de Sitter case, the model still under- 
predicts the faint end of the source count slope in the op- 
tical, UV and near-IR, though to a lesser degree. There is 
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5.4,n = 1.3, A 



Vo 



0.4 [solid curve]) and semi- 



analytic differential luminosity functions (Kaviani et al., 2002 
[dashed curve]) at redshifts, z=0.0, 1.0 (offset +1), 2.0 (offset 
-1-2.5), 4.0 (offset 4-4.5). The two models are broadly consistent. 



also an over-prediction of very faint submillimetre counts 
and the high-redshift SFR. These effects are not removed 
by accounting for ellipticals. It is possible that further cor- 
rections to the model could resolve these problems. The A 
cosmology does not suffer from these effects, though it does 
over-predict the FIR region of the background spectrum. 

(iv) The mean local value of Ay of the best-fitting model 
is 0.4, which is consistent with what we know about the 
opacity of the local universe. 

(v) The form of the luminosity function and the luminos- 
ity and density evolutions is broadly consistent with those 
predicted by a semi-analytic model. 

(vi) The predicted form of the SFH provides approxi- 
mately the same goodness of fit in the E-dS case, and a 
better fit in the A case. 

(vii) It is not possible to to reproduce the multi- 
wavelength counts using a model without density evolution 
(the monolithic model). 

(viii) The model provides a fairly good fit to data from 
UV to submillimetre. Future work will investigate expanding 
this to the radio and X-ray regions. 
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